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In this paper we consider an exactly solvable model which displays glassy behavior at zero 
temperature due to entropic barriers. The new ingredient of the model is the existence of different 
energy scales or modes associated to different relaxational time-scales. Low-temperature relaxation 
takes place by partial equilibration of successive lower energy modes. An adiabatic scaling solution, 
defined in terms of a threshold energy scale e* , is proposed. For such a solution, modes with energy 
e 2> e* are equilibrated at the bath temperature, modes with e <C e* remain out of equilibrium 
and relaxation occurs in the neighborhood of the threshold e ~ e*. The model is presented as a 
toy example to investigate conditions related to the existence of an effective temperature in glassy 
systems and its possible dependence on the energy sector probed by the corresponding observable. 



I. INTRODUCTION 

The study of exactly solvable models has been always an active area of research in the field of statistical physics. 
They help us to grasp general principles governing the physical behavior of realistic systems which, due to the 
complicated interactions among the different constituents, cannot be predicted using standard perturbative techniques. 
Glasses in general are systems falling into this category. The slow relaxation of glasses observed in the laboratory is 
a consequence of the simultaneous interplay of its constituents that yields a very complex and rich phenomenology. 

It is well known that glasses fall out of equilibrium when the characteristic observation time is larger than their 
relaxation time. Because the relaxation time is strongly dependent on temperature, it turns out that glasses are 
immediately out of equilibrium as soon as the temperature is few degrees below the glass transition. Well below the 
glass transition temperature T g no time evolution is apparently observed in the glass and one is tempted to conclude 
that the glass is in a stationary state. Nothing more far from the truth. Glasses still relax but slow enough for 
any change to be observable in laboratory time-scales. Old experiments on polymers reveal that the slowly relaxing 
state corresponds to an aging statetl. That is, if the system is perturbed while being in its aging state, then the 
characteristic time associated to the response of the system scales with the age of the system (i.e. the time elapsed 
since it was quenched) . Another way to look at this aging phenomena is to evaluatei-the time autocorrelation function. 
It is observed that the typical decorrelation time scales with the age of the systemo. 

A simple scenario to explain these results is the following. Consider a liquid well above T g where correlations decay 
exponentially with time. One may consider the resultant behavior of the liquid as the superposition of different and 
independent harmonic modes. Each of this energy modes corresponds to a normal mode of a system and describes a 
collective oscillation of N atoms around their local minimum. This is the harmonic approximation which is known 
to work quite well in liquids. Nevertheless, already as the temperature goes below a critical temperature T c (the 
transition temperature predicted by the mode-coupling theoryB) other collective modes different from the standard 
vibrational ones become important. The nature of these modes is quite different from the usual harmonic normal modes 
because they do not represent oscillations around a given configuration within a metastable well but transitions among 
different wells. These modes are reminiscent of some type of instanton solutions recently computed in the framework 
of some spin glass modelsu. Below T c relaxational dynamics proceeds by activation over the barriers characterizing 
these collective modes. Now, the main difference between these collective modes and the usual harmonic modes 
relies on how they relax to equilibrium when put in contact with a thermal bath at temperature T. Relaxation to 
equilibrium is determined by the height of the energy barriers separating different modes. Suppose a given normal 
mode has frequency u>k and energy oc uj\. The relaxation time for each of these modes is typically of order 
Tfc ~ exp(i?/ c / KbT). Therefore, as the energy is lower than the thermal bath temperature this mode rapidly 
equilibrates. On the contrary, if Ek 3> KbT, this mode remains frozen. Collective modes arc different. As the 
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reference energy of the collective modes deplete, the typical barrier separating these modes increases leading to the 
contrary behaviour and to super-activation effects. While high energy collective modes are separated by low barriers, 
low energy collective modes are separated by high barriers. A simple schematic representation of this scenario in a one 
dimensional configurational space is shown in figure 0. This behavior is common to the majority of exactly solvable 
glassy modelg3~LJ , phenomenological trap modelsa and kinetically constrained modelsa to cite a few. 




{configuration space} 



FIG. 1. A one dimensional example of an energy landscape of collective modes: high energy collective modes are separated 
by low barriers while low energy collective modes are separated by high barriers. 

In what follows we will use the generic word mode to refer to this kind of collective excitation. Let us label the 
modes with the integer variable r and let us denote their energy by e r . Let us suppose that the energy levels are 
ordered from lower to higher energies according to the label r. It is natural to assume that there is a characteristic 
mode r* with associated energy e* such that, all modes with e r e* have already relaxed, while in the other limit, 
e r •C e*, all modes are frozen. If the system is quenched well below T c then equilibration cannot be achieved in 
laboratory time scales, this means that all modes below e* remain frozen while modes above e* remain equilibrated 
at the bath temperature. The energy threshold e* decays with time because, as time goes by, higher barriers are 
accessible to the system. The resulting scenario is that of a liquid where collective modes above e* are in some sort 
of local equilibrium at the temperature of the bath T while modes below e* are frozen. This scenario, as it stands, 
is too naive because it is based on the assumption that there are no dynamical correlations between the different 
modes, i.e. high energy modes do not influence low energy modes. While this is true in the equilibrium state it may 
not be valid (an indeed it is not) when any type of local dynamics induces correlations between different modes. Still 
it is interesting to investigate in which conditions. such a description turns out to be correct. Examples where these 
type of description holds are mean-field modelsEl'E£rtj. By mean-field we mean those models where there is no spatial 
dimensionality associated with the set of interactions. In these cases, different modes are not related to different 
length scales. Therefore local dynamical rules do not necessarily induce correlations between the different modes. 
Generally speaking, the identification of this energy scale remains an open problem for which we do not have yet a 
complete understanding. 

This energy scale e* is related to what_has. .received the name of effective or Active temperature T c g in the most 
recent literature about glassy dynamicstj EJEjrj. T G ff is an effective time-dependent parameter describing equilibrium 
fluctuations for those thermalized modes with e ^ e* = ksTcs. While a thermometer coupled to the fast modes is 
expected to measure the bath temperature, a thermometer coupled to the slower ones should yield the effective 
temperature. Hence, the effective temperature corresponds to an energy threshold which separates collective modes 
which are frozen from those which have relaxed to the bath temperature. Relaxation, then, takes place at energy 
scales in the neighborhood of the threshold value e* . For real systems this energy threshold is related to the typical 
volume of the drop which is able to release strain energy during its relaxation to equilibrium. A scenario on this type 
of physical mechanism for glassy relaxation has been recently introduced inEd where this threshold energy e* has been 
related to the size of the cooperative region as an explanation of the fragility and super-activation anomalies in real 
glasses. 

The existence on an effective temperature is tightly related to the validity of some approximations used in the context 
of slowly relaxing systems such as the adiabatic approximation. Within the adiabatic approximation one obtains 
a Markovian description for the dynamics. Again, this Markovian description encodes within a single parameter 
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(the effective temperature) all the complicated previous past of the system. The adiabatic approximation has been 
shown-to-give the correct asymptotic dynamical behavior in some simple models of glasses such as the Backgammon 
modeE3'B2l while for the majority of the other most famous models in the literature (for instance the p-spin models) the 
implementation of such approximation is generally not yet known.— Quite probably a Markovian description for glassy 
dynamics is unrealistic and the original idea of experimentalistsc.3 to encode the dynamical behavior into effective 
parameters such as the effective pressure or the effective (else said fictive) temperature could be only labels without 
a deep physical meaning. 

These ideas have been contrasted in several exactly solvable models. All these models have the advantage of 
being mean-field, hence dynamical equations can be closed in one or another way. Still, the mean-field character 
of these models does not generally allow to investigate in a simple way the relaxational properties of the different 
modes of the spectrum. In general, the dynamics of all these models has been studied by addressing the closure 
of the dynamics associated to global quantities (such as the energy or correlation-response functions) which do 
not discern the G^ffitribation to the global Kelp;atipii of the different energy modes. Examples of these mod&ls-s*e 
the p-spin modeltHi3li3't3, the BackgammonE3£3SEa the harmonic oscillator and other spherical spin models&E^LZl. 
The question of the existence of an effective temper ature i s also tightly related to the particular way in which 
the fluctuation-dissipation theorem (FDT) is violatedl3E3EZl. The question whether there exists a single effective 
temperature describing the violation of FDT is still controversial (see, for instance,C3) . We believe that this and other 
related questions can be better addressed, even at the level of mean-field models, by analyzing the contribution of 
all different modes. The analysis of the harmonic vibrational modes in glasses is already an interesting step in this 
directions although here we have in mind other type of collective modes. 

The purpose of this paper is the study of a new model where the dynamical relaxation of the different energy modes 
can be made explicitly clear. This model is what we refer to as the disordered Backgammon model (DBG model) and 
consists in a generalization of the Backgammon model to allow different energies for different boxes. Again, just like 
its predecessor, the slow relaxation of this model is due to entropic barriers. For the DBG model we show the existence 
of an energy threshold e*(t) which separates equilibrated form non-equilibrated modes. The DBG model (with its 
associated threshold e*(t)) provides a microscopic realization that is reminiscent of some rphenomenological models 
proposed in the past such as the trap model in a tree considered by Bouchaud and DeanEJ. The advantage in the 
DBG model is that now one can exhaustively investigate the distinct relaxation of each of the different energy modes 
verifying whether the scenario of the effective temperature presented before holds. This will help to better understand 
the meaning of the effective temperature and its relation with the violation fluctuation-dissipation theorem. 

In section II we introduce the disordered Backgammon (DBG) model. Its thermodynamic properties are reported 
in section III, where the special case of very low temperatures is explicitly worked out. In section IV we present the 
dynamical equations whose solutions are found within a new type of adiabatic approximation in section V. There the 
features of such an approximation are carefully analyzed. In section VI numerical results are presented for two specific 
models belonging to the family of DBG models . Finally, in section VII, a novel method is introduced in order to 
estimate the threshold energy scale e* directly from the dynamics. In section VIII we present our conclusions. Some 
technical issues are presented in four different appendices. 

II. THE DISORDERED BACKGAMMON MODEL 
A. Definition of the model 

Let us take N particles which can occupy N boxes, each one labeled by an index r which runs from 1 to N. Suppose 
now that all particles are distributed among the boxes. A given box r contributes to the Hamiltonian with an energy 
—e r only when it is empty. In this case the total Hamiltonian of the system reads 

N 

Ti = -^er<5n r ,o , (1) 

r=l 

where S is the Kronecker delta and n r denotes the occupancy or number of particles in box r. The e r are quenched 
random variables extracted from a distribution g(e) which we assume to be defined only for e > 0. The interest 
of this definition will be discussed below. Like in the original Backgammon (BG) modeB we consider Monte Carlo 
mean-field dynamics where a particle is randomly chosen in a departure box and a move to an arrival box a is 
proposed. If d denotes the departure box, the proposed change is accepted according to the Metropolis rule with 
probability W(AE) = Min(l, exp(— (3AE)) where AE = e a &n a fl — £d$n d ,i- Note that the departure box satisfies 
rid > 1 and departure boxes are chosen with probability n^/N. This dynamics corresponds to Maxwell statistics 
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where particle are distinguishable and differs from the one corresponding to Bose statisticfH ! c3E3£3 where particles 
are indistinguishable and arrival boxes are chosen with uniform probability l/N. 

In the dynamics the total number of particles is conserved so that the occupancies satisfy the closure condition 



Now we want to show that the interesting case corresponds to the situation where g(e) is only defined for e > 0. 
In this case, the dynamics turns out to be extremely slow at low temperatures, similarly to what happens for the 
original BG model. The difference lies in the type of ground state. The ground state of (|l|) corresponds to the case 
where all particles occupy a single box, the one with the smallest value of e. Let us denote by eo this smallest value. 
Then the ground state energy is given by the relation Eqs = — J2r=i e »' + e o- Since all the e are positive no other 
configuration can have a lower energy. If 17(e) is a continuous distribution the ground state is also unique. Now it is 
easy to understand that, during the dynamical evolution at zero temperature, all boxes with high values of e become 
empty quite soon and the dynamics involves boxes with progressively lower values of e. The asymptotic dynamics is 
then determined by the behavior of the distribution g(e) in the limit e — > 0. If <?(e) ~ e Q , for e — > 0, we will show 
that the asymptotic long-time properties only depend on a. Note that the normalization of the g(e) imposes a > — 1. 
This classification includes also the original BG model where there is no disorder at all. In that case g(e) = 6(e — 1) 
so the distribution has a finite gap at e = 0. The behavior corresponding to this singular energy distribution can be 
obtained from the previous one in the limiting case a — > 00. 

One important aspect of the model (Q) is that, in the presence of disorder, it is not invariant under an arbitrary 
constant shift of the energy levels. Actually, by changing e r — > e' r = e r + c with c > 0, the model turns out to be 
a combination of the original model plus the classical BG model. After shifting, the new distribution g(e' — c) has 
a finite gap (equal to c plus the gap of the original distribution) . The new model corresponds again to the a — > 00 
case and the asymptotic dynamical behavior coincides with that of the standard BG model. As we will see later the 
present model is characterized by an energy threshold e* which drives relaxation to the stationary state. Only when 
the energy threshold can go to zero we have a different asymptotic behavior. For all models with a finite gap, e* 
cannot be smaller than the gap, hence asymptotically sticks to the gap and the relaxational behavior of the DBG 
model with a finite gap corresponds to that of the standard BG model. 

One of the outstanding features of this model is that a description of the dynamics in the framework of an adiabatic 
approximation turns out to be totally independent from the type of distribution g(e) (and hence on a) despite of the 
fact that the asymptotic long-time behavior of the effective temperature and of the internal energy depend on the 
value of a. 



N 




(2) 



B. Observables 



Like in the original BG model we define the occupation probabilities, Pk, that a box contains fc particle 




(3) 



and the corresponding densities that a box of energy e contains fc particles 




fc > ; 



(4) 




(5) 



The Pk and the gu are related by 




(6) 



and the conservation of particles reads 
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fc=0 k=0 

The energy can be expressed in terms of the density go(e) as 



OO OO , \ 

£> = 1 ; X>(e) = r J [£ ,\, =9(e). (7) 
— ' Jo 5(e)* 



E=- dee g (e) . (8) 
Jo 

This set of observables depend on time through the time evolution of the occupancies n r of all boxes. In the next 
section we analyze the main equilibrium properties of the model. 

III. EQUILIBRIUM BEHAVIOR 

The solution of the thermodynamics proceeds similarly as for the case of the original BG model. The partition 
function can be computed in the grand partition ensemble. It reads 

OO 

Z GC = Z c (N)z N , (9) 

N=0 

where z = exp(Pn) is the fugacity, fx is the chemical potential and Zq{N) stands for the canonical partition function 
of a system with TV particles. The canonical partition function can be written as 

N ATI N 
Zc= Y.TTn~ — : ex P(P'52 e r S T+fl)8 N y „ r . (10) 

n r =0 llr=l " ' 



r=l 



where Sij is the Kronecker delta. Introducing this expression in (^|) we can write down Zqc as an unrestricted sum 
for all the occupancies n r 



oo N 

Z 



Z GC = N\ II — f exp(/3e r 5 nr)0 ) . (11) 



nJ- 

n r —0 r—1 



The factor N\/Y\ r n r \ in the partition function Zq is introduced to account for the distinguishability of particles. 
This factor leads to an over-extensive entropy (i.e., the Gibbs paradox) which can be cured eliminating from Zqc the 
over-counting term AH in the numerator. The final result is 



f N 

Z GC = exp < 2J log 

I r=l 

yielding the grand-canonical potential energy per box 

JV 



z 

E— exp (/3e r <5„ !0 ) 



n=0 



(12) 



G=F-n=-T 



log(Zac) 



N 

r=l 



T V log [exp(/3e r ) + exp(z) - 1] = -T / g(e) log [exp(/3e) + exp(z) - 1] de , (13) 



where F is the Helmholtz free-energy per box. The fugacity z is determined by the conservation condition (g) which 
reads 

dG dF 

— — = — 1 or, equivalently, — — = , (14) 

d[l O/J, 

yielding the closure condition 

9(e) 



, \ 7* — TT ■ (15) 

o exp(pe) + exp(z) — 1 zexp(z) 

This equation gives the fugacity z as function of (3 and, from (|l^) and its derivatives, the whole thermodynamics. 
In particular, the equilibrium expressions for g%(e) are 
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eg/ \ = z fc g(e)exp(/3e4,o) , , 

5fe W fc!(exp(/3e) + exp(z) - 1) ' 1 J 

the corresponding being given by eq.(^J) which, together with the closure relation (|l5|), leads to the expression 

= 6 k , (l ~ eMZ \~ 1 ) + (1 - S^j^—n ■ (17) 
.zexp(z) fc!exp(z) 

Starting from eq.(|[) the equilibrium energy density is obtained as 

E°« = -r d e l 9 [ £)eM ^ : ■ (18) 
7 exp(/3e) + exp(z) - 1 

All these expressions can be evaluated at finite temperature. Note that, although the values of P^ q in ( [i"7| ) are 
independent of the disorder distribution g(e), they directly depend on that distribution through the equilibrium value 
of z (which obviously depends on the 17(e)). Of particular interest for the dynamical behavior of the model are the 
low-temperature properties that we analyze below. 



A. Thermodynamics at low temperatures 

A perturbative expansion can be carried out close to T — > to find the leading behavior of different thermodynamic 
quantities. Let us start analyzing the closure condition (|l5|). Doing the transformation e' = [3e, eq. (|l5|) can be 
rewritten as 

r r«_ t . ' . (19 , 

Jo exp(e') + expyz) — 1 zexp(z) 

In the limit T — > the fugacity z depends on the behavior of g(Te') in the limit T = 0, i.e. on the behavior of g(e) 
for e — * 0. Assuming g(e) ~ e a for e <C 1 we define the function h(e) through the relation g(e) = e a h(e), where h(e) is 
a smooth function of e with h(0) finite. The integral can be expanded around T = by taking successive derivatives 
of the function h: 



J expe' + exp z -If- k\ 

K — \J 



Using the asymptotic result z — > 00, when T — * 0, everything reduces to estimate the following integral in the large z 
limit: 

dx / n XC ' +k , \ 7 ~ z ° +k+1 exp(-z) . (21) 

exp(x) + exp(z) - 1 y J y J 

The term k = in the series yields the leading behavior for z, which turns out to be 

z ~ /3§T3 (22) 

In a similar way the energy can be computed to leading order in T, 

E = E GS + aT + 0(T 2 ) , (23) 

giving a finite specific heat at low temperatures. In section VI we will show explicitly such a behaviour for two specific 
DBG models, one with a = 1 (defined in eq. (^|), sec. VI) and the other with a — (eq. (fj9|) , sec. VI). Solving 
eq. ( |l5| ) numerically for z(T) in each specific model and inserting z(T) in the expression (|18D for the equilibrium 
energy density, we get the energy dependence on the temperature. In all cases, at equilibrium, the energy is linear 
for low T, as predicted in eq. (j23|). It yields, therefore, a finite specific heat as corresponds to a classical model with 
Maxwell-Boltzmann statistics (see fig. (|])). 
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FIG. 2. Left: energy dependence on temperature at equilibrium for two DBG models, with a = 1 and a — 0. Also the 
energy for the standard BG model (a — > oo) is plotted. Right: specific heat for the same models. For all cases (disordered or 
not) it turns out to be finite at T — 0. 



IV. DYNAMICAL EQUATIONS 



Here we consider the dynamical equations for the occupation probabilities Pk and their associated densities <7fc(e). 
The dynamical equations in this model are derived in a similar way as for the standard BG model. The main difference 
is that in the DBG the equations for the occupancies probabilities Pk do not generate a closed hierarchy of equations. 
Only for T = such a closed hierarchy is obtained. As we will see later, this has important consequences when we 
become interested in the zero-temperature relaxation. 

A hierarchy of equations can only be obtained at the level of the occupation probability densities gfc(e). A detailed 
derivation of these equations is reported in the Appendix A. Here we show the final result, 



dgo(<0 
dt 



dgi(e) 
dt 



ggfc(f) 

dt 



1 



de 



'go(e') (e 



1 



* + Pl (1 - e -/3e) + f de' gi (e') (e~^-^ - l) 



2 52 (e)ll+/ deg (e)e-^-P Q \ - g x {e) 



de'go(e') (V^'-> - 1 



5o (e) 



+ Pi (1 - e-^) + de'gi(e') (e~^-^ - l) 



(k + l)flk+i(e) [ 1 + jf deg (e)e-^ - P a 
1 /,' + /■■( J d'U-.Acu " /' 



■fffc-i(e) 



fc = 1 



fc > 1 



(24) 



(25) 



(26) 



The equations for the Pk are directly obtained by integrating the gfc(e) according to (|^). The result is 
dP k {t) 



dt 



(k + 1) {P k +i(t) - P k (t)) + Pk-i + (J deg (e)e-^ - P ) - 4,o - kP k (t) + (k + l)P k +i(t)] (27) 



with P_i =0. It is easy to check that the equilibrium solutions (^6() are indeed stationary solutions. As previously 
said, for general /3, the equations for the Pk do not generate a hierarchy by themselves but depend on the <7fc(e) through 
the distribution go(e) in (|2^). Nevertheless, a remarkable aspect is that they generate a well defined hierarchy at T = 
which coincides with the equations of the original BG modelEj. These are 



dP k {t) 
dt 



= {k + 1) {P k +i(t) - Pk{t)) + Pk-i - Pa [4,1 - 4,0 - kP k {t) + (k + l)JVu(t)] 



(28) 



7 



It is easy to understand why at T = the dynamical equations are independent on the density of states g(e). The 
reasoning is as follows. For T = 0, all moves of particles between departure and arrival boxes with different energies 
€d and e a depend on the precise values of these energies only when the departure box contains a single particle and 
the arrival box is empty. But such a move does not lead to any change in any of the Pk , hence dynamical equations 
for the Pk remain independent of g(e). Obviously this does not hold for other observables such as the energy E (i.e. 
the mean value of e over the density go(e)) and higher moments of 3fe(e). 

This observation is crucial, since we look at the glassy behaviour at T = 0. It turns out that the analysis of the 
dynamical equations for the DBG at T = decomposes into two parts. On the one hand, the equations for the Pk 
coincide with those of the original non-disordered BG model. On the other hand, in order to analyze the behavior of 
the energy one must analyze the behavior of the hierarchy of equations for the <7fe(e) which is quite complicated. In 
the next section we will see how the analysis of these equations can be done within the framework of a generalized 
adiabatic approximation. 



V. THE ADIABATIC APPROXIMATION 



A. Standard adiabatic solution for the Pk 



In this section we are interested in the solution of the dynamical equations at T = 0. We already saw that the 
equations (^8|) for the Pk coincide with those of the original non-disordered BG model. Consequently, the same 
adiabatkj-approximation used for the Pk in the original BG model is still valid for the DBG. Let us remind the main 
resultsE3'E2l. The key idea behind that approximation is that, while Pq constitutes a slow mode, the other Pk with 
k > are fast modes. Hence they can be considered as if they were in equilibrium at the hyper-surface in phase space 
Pq = constant, this constant being given by the actual value of Pq at time t. In the original BG model Pq = —E, 
hence thermalization of the fast modes Pk (k > 0) occurs on the hyper-surface of constant energy. For the DBG 
model this is not true, the hyper-surface where equilibration of fast modes occurs does not coincide with the constant 
energy hyper-surface simply because the energy and Pq arc different quantities. Indeed, we will see later that their 
leading time behavior is different . 

At T — the equation for P (eq. (||) for k = 0) reads 

BP 

-^=P 1 (1-P ). (29) 

If local equilibrium is reached on the hyper-surface of constant Pq we can relate Pi to Pq using eqs.(^). The 
simplest way of dealing with ( p9| ) is to relate both Pi and Pq to the time-dependent fugacity z* writing down a 
dynamical equation for z*. Using eq.(|l7|) this yields 

dz* = z*(exp(z*) - 1) 

dt exp(z*)(exp(z*) - z* - 1) ' 1 ' 

In the large time limit z* diverges and the leading asymptotic behavior is given by z* ~ log(t) + log(log(f)). The 
occupation probabilities are, then, given by the relations (O) replacing z by the time-dependent fugacity z*, as long 

as z* > 1, 

ft - (l - i) + (1 - tffc .o) ^f^J • (3D 
Hence, according to (p^), the inverse effective temperature is given by the relation 

fcg ~ (z*)%& ~ (log(t))^ (32) 

and the effective temperature depends on the properties of the disorder distribution g(e) in the limit e — > through 
the value of the exponent a. Clearly, when the density of levels decreases as we approach e = 0, the relaxation turns 
out to be slower; the limiting case being the original BG model for which a —> oo and /3 e a ~ log(f). In the other limit 
a — > — 1, when disorder becomes unnormalized, the inverse effective temperature diverges very fast. Already from 
( p3| ) one can anticipate that the same asymptotic behavior holds for the energy (see eq. (jij), hence a interpolates 
between fast relaxation (a = —1) and very slow relaxation (a = oo). A relaxation slower than logarithmic is not 
possible in the present model. 
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B. Generalized adiabatic solution for the <?&(e) 



The equations for the gfc(e) at T — are 



dffo(e) 

at 

3gi(<0 
dt 

dgkje) 

dt 



9i(e) 



de'g (e') 



g (e) J de' gi (e') fe = , (33) 

de'g (e') + g (e) / de' 9l (e') k = 1 , (34) 



2< ?2 (e)(l-P )-< ?1 (e) 

{k + l)g k+1 (e){l-P )-g k (e)[l + k(l-P )}+g k . 1 (e) k>l. (35) 



To solve the dynamical equations for the g k {^) in the adiabatic approximation we note that, contrarily to the 
global quantities P k , they cannot be equilibrated among all different modes. The reason is that, due to the entropic 
character of the relaxation, very low energy modes are rarely involved because the time needed to empty one further 
box increases progressively as time goes by, hence they cannot be thought as effectively thermalized. Note that in the 
original BG model all boxes have the same energy, hence there is a unique class of modes. For the general disordered 
model we expect the existence of a time dependent energy scale e* separating equilibrated and non-equilibrated 
modes. The mechanism of relaxation is the one we reminded in the introduction when speaking about the behavior 
of collective modes. At zero temperature there is no thermal activation and the equilibrated modes are in the sector 
e»6* while the non-equilibrated modes are in the other sector e e* . The value of e* can be easily guessed. After 
quenching to zero temperature the system starts to relax to its ground state. Because there is no thermal activation, 
relaxation is driven by entropic barriers, i.e. flat directions in configurational space through which the system diffuses. 
Entropic relaxation is energy costless so its rate is determined by the number of available configurations with energy 
smaller or equal to the actual energy. A simple microcanonical argument giving the relaxation rate goes as follows. 
Let us denote M the number of occupied boxes and f2(M) the number of configurations with M occupied boxes. The 
typical time to increase by one unity the number of empty boxes is given by, 

-~ ^ (36) 



Q(M - f ) 

For distinguishable particles we have f2(M) ~ M N . In the large N,M limit we get 



/ M ,n ,N. 

fc) ^m } - (37) 



exp -- , (38) 



Using the relation Pq = 1 — M /N and (^) we find 

dP AM _ 1 
~~dT ~ ~ NAt ~ t I -Pi 

where we have used the fact that at zero temperature the number of occupied boxes can only decrease by one unity, 
thus AM = -I. 

This yields the result Pq ~ f — i og ( t ) + iog(iog(t)) m agreement with the adiabatic approximation. Actually, from the 
solution in (^o|) for z* and using the adiabatic relation (|3f"|), we obtain Pq ~ f — 1/z*, which yields the same result. 



The typical relaxation time (37) behaves, then, like r ~ exp( 1 _ P(j ) — exp(z*). 

If the threshold e* plays the role of an energy barrier and f3 g accounts for the effective thermal activation due to 
entropic effects, we obtain, for the typical relaxation time, r ~ exp(/3 c ffe*). This expression is only valid to leading 
order. As we will see below there are sub-leading corrections to this expression arising for the fact that the relaxation 
time is better described by the expression r ~ exp(/3 e ffe*)/(/3 c ge*) (see eq. (|45|)). Hence, at a given time-scale t (i.e., 
the time elapsed since the system was quenched) all modes where rCt are equilibrated at zero temperature (which 
in this case is the temperature of the thermal bath) and therefore frozen. Modes with T>i, although dynamically 
evolving, are also frozen because the barriers (in this case entropic barriers) are too high to allow for relaxation within 
the time-scale t. 

Only those modes whose characteristic time is t ~ t are relaxing at a given time-scale t. We get for the time 
dependent energy scale e* and the effective temperature the relation e* ~ log(t)//3 c g and this yields the leading 
behavior 

e* ~ (log(t))-=k . (39) 







According to what has been said one can impose the following Ansatz solution for the g k {e)- If g cq stands for the 
equilibrium density at T = (i.e., according to eq. (|l6|), g^ q = g{e)Skfi) then we have, 



Ag k (e) = g k (e) - g^(e) = ^g k & (40) 



e 



where APk = Pk ~ P k q — Pk — Sk.Q and gk(x) decays pretty fast to zero for x > 1. This expression tells us the 
following. Above e* the gk{e) have relaxed to their corresponding equilibrium distributions at the temperature of the 
bath (in this case the bath is at zero temperature). On the other hand, in the sector of the energy spectrum where 
e < e* , the densities gk are still relaxing (specially in the region e/e* ~ 0(1)). Since the relaxation is driven by the 
shift in time of the threshold energy e* the proposed scaling solution Ansatz seems quite reasonable. The prefactor 
APk /e* is introduced to fulfill condition (||) . Furthermore the condition dx g k (x) = 1 is imposed on the scaling 
function <jk- 

In the Appendix B we show how this Ansatz closes the set of equations ( |33| ) reproducing also the leading asymptotic 
behavior for e* and z* which turns out to be, 

' (41) 



(log(i)) — 

l^logW (42) 



For later use we define the following function 

, W _£gjf£_ fc £), (43) 

which scales as function of e/e*. The scaling relation ( flof ) yields the leading asymptotic behavior of all observables 
different from the occupation probabilities Pk- For instance, the energy is given by E = — J Q de e go(e); using the 
scaling relation ( pio| ) and the asymptotic expression ( f4l| ) we get for the leading term 

r°° i 

E-E GS ~- de [ 50 (e) - g(e)} e ~ (e*) a+2 — . (44) 



log(i)- 



+ 2 
+1 



Note that the asymptotic scaling behavior of the energy is the same as for the effective temperature T c g = 1 //3 e ff 
in agreement with the quasi-equilibrium hypothesis (see (|23|)). An important result is that the threshold e* decays 
slower to zero than the effective temperature. A case where this difference can be clearly appreciated corresponds to 
the case where the density of states vanishes exponentially fast g(e) ~ exp(— A/e). In this case e* decays slower than 
logarithmically, namely like 1/ log(log(i)) (see appendix B for details). 



C. Relaxational spectrum in equilibrium 



One of the crucial features behind the applicability of the adiabatic approximation is that the long-time behavior at 
zero temperature finds its correspondence with the low-temperature relaxational properties of the equilibrium state. 

To analyze the spectrum of relaxation times T eq (e) in equilibrium we expand up to first order in perturbation 
theory the dynamical equations for the gk(e) around their equilibrium solutions g c k q ( e )- Using the expansion gk{e) = 
<7^ q (e) + 8gk(e) we get a set of equations for the variations Sgk{e). These are shown in the Appendix C. 

A complete derivation of the relaxation time r(e) in equilibrium is complicated. But it is easy to convince oneself 
that the relaxation time is asymptotically (in the limit T — > 0) strongly peaked around the threshold energy e*. For 
e 3> e* the relaxation time is small because the population of high energy boxes in equilibrium is rather small. On 
the other hand, for e/e* <C 1 the relaxation is estimated to be finite and independent of T. This result is derived in 
the aforementioned Appendix C where we show that the maximum relaxation time occurs for e around e* . Starting 
from eqs. ( |6^ ) for Sgo(e) and Sg\(e) and making use of the adiabatic Ansatz (pfl|), we find, for e ~ e* , 

r(e*) ^ ^ (45) 

where e*(T) ~ T 1 ^ 2+a \ is the asymptotic temperature dependence of the threshold energy at low temperature. This 
yields for the temperature dependence of the relaxation time, 
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Q + l 

13 &n 

r(T) ~ (46) 
/?°+ 2 

showing there is activated behavior as a function of the temperature but with a relaxation time which increases 
slower than Arrhenius as T 0. Note that for the standard BG model corresponding to a — > oo we obtain an 
Arrhenius behavior and in the opposite limit a — > — 1 the relaxation time does not diverge anymore. 



VI. NUMERICAL RESULTS 



In this section we numerically check the main results obtained in the previous sections. In particular, we want 
to show the existence of the threshold energy e* separating equilibrated from non-equilibrated energy modes. We 
have compared three different models characterized by three different types of distributions (figure (|3|)). All three 
distributions were chosen to satisfy the conditions 

I* oo />oo 

/ de g(e) = dee g{e) = 1 , (47) 
Jo Jo 

in such a way that the ground state has energy Eqs = — 1 in the limit N — > oo for all three cases. The models are 
the following ones: 

• Case A: Non disordered model with a gap (figure (|a)). This is the original BG model where g(e) = S(e — 1). 
This case corresponds to a — > oo, therefore e* = 1 and the threshold energy is time independent. The energy is 
expected to decay like E + 1 ~ T e g ~ l/log(t). As previously said in section II, the same behavior is expected 
for any disorder distribution g(e) with a finite gap. 

• Case B: Disordered Model without gap but g(0) — (figure (||b)). We have considered the distribution 

ff ( £ ) = | e exp(-J e 2 ) . (48) 



This case corresponds to a = 1. The energy threshold e* scales like l/-\/log(i) and the effective temperature 
and the energy scale like E + 1 ~ T c s ~ l/(log(£))i. 

Case C: Disordered Model without gap and g(0) finite (figure (j^c)) . We have considered the distribution 

3 (e) = -exp(--) . (49) 

This case corresponds to a = 0. The energy threshold e* scales like 1/ \og(t) and the effective temperature and 
the energy scale like E + 1 - T cS - l/(log(t)) 2 . 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 



e e e 

FIG. 3. Probability distribution of the energy weights of the boxes of the DBG model, (a) The standard Backgammon model 
has no disordered distribution, all boxes have the same weight, (b) The probability distribution function of a DBG with a = 1, 
at very low energy the density of boxes goes to zero, (c) A second DBG model with a = 0. Here the probability of having 
boxes with energies arbitrarily close to zero is finite. 



11 



In figure (||) we plot the decay of the energy for all three models. Simulations were done for N = 10 4 , 10 5 , 10 6 
boxes (the number of particles is identical to the number of boxes) showing that finite-size effects are not big in the 
asymptotic regime. We show data for one sample and N = 10 6 . We plot the energy as function of time starting from 
a random initial condition (particles randomly distributed among boxes, E(t = 0) = — 1/e). As clearly seen from the 
figure relaxation is faster for the Case C and slower for the standard BG model (Case A). 




t 

FIG. 4. Energy as a function of time for the three models discussed in the text. The energy was computed averaging over 
time intervals At = 2 k , k being an integer number. The lower curve represent the relaxation in the a — DBG model (C), the 
middle curve the a — 1 DBG model (B) and th upper curve the standard BG model (A). 

The different asymptotic behaviors are shown in figure (g). There we plot (E(t) — Eas)(\og(t)) x with A = §if • To 

avoid finite-size corrections when the energy is close to its ground state we computed exactly Eqs — w(~ Sr=i £r " 4 ~ 
£min) where e m i n is the minimum value among the e's. The different curves saturate at a finite quantity corresponding 
to the asymptotic leading constant. Note that convergence is slow, showing the presence of subleading logarithmic 
corrections to the leading behavior. 



2.5 




t 

FIG. 5. (E — _Ecs)(log(t)) A plotted as function of time with A = 2±2 for the three different models discussed in the text. 
The upper curve refers to case C, the middle one to case B and the lower one to standard BG model (case A). 

Let us now analyze the shape of the probability densities gu (e) ■ For these distribution we also only show results 
for N = 10 6 because a smaller number of boxes results in higher noise in the curves. The distribution probabilities 
were numerically computed by binning the e axis from e = up to e = e max where e rnax is the maximum value of 
e r among all the N boxes. One hundred bins are enough to see the behavior of the time evolution of the different 
distributions. In figures (||) and (]?]) we show the <?o( e ) f° r cases B and C respectively. Note that the go(e) converge to 
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the asymptotic result g(e) for e > e* in agreement with the adiabatic solution (40) while they are clearly different for 
e < e*. The value of e* where go(e) deviates from the asymptotic curve g(e) shifts slowly to zero (like l/(log(i))2 or 
l/log(i) for cases B and C respectively), as can be seen in figures (§J7|)). 




c 
so 




2 3 
e 

FIG. 7. Distribution g (e) for case C for different times 2 fe with k = 4, 6, 8, 10, 12, 14, 16, 18, 20 (from bottom to top). The 
continuous lines is g(e) given in eq. Q49|) . 



We do not show results for the other <?& (for instance g±) because they decay very fast to zero (already for t = 2 17 
there are no occupied boxes with more than one particle). Instead, in figures (^) and @ we verify the adiabatic 
Ansatz, eqs.( ^0| , ^3| ), for the densities go and gi in the two models B and C. Figure (||) plots Go(e) for both models. 
Figure (^|) plots Gi(e) for both models. We have used the relation ([39|) together with z* = log(t) + log(log(t)) yielding 
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n , v A g0 e (log t + log log t ) e 

G (e) = = go(-) (50) 

log(t) "+ 1 e 

Gl(e) = ^i(^ = ^(±) (51) 
U 7 log(i)— * V e*> V y 

Note that the scaling is pretty well satisfied and that the gk{x) indeed vanishes for x ~ 1 yielding an estimate for 
e* is both cases. We find, e* ~ 6/ ^/log(t) for case B and e* ~ 12/log(t) for case C. Note also that the quality of 
the collapse of the Go is slightly worse for case B than for case C. We think this is due to the stronger subleading 
corrections to the shift of e* which decays slower to zero for case B. Hence the asymptotic regime is reached only 
for later times. Indeed, as figure (^|) shows, the value of e* obtained within our time-scales has not yet reached 
the maximum of the distribution g(e), so that we are still far from the asymptotic behavior g(e*) ~ e*. Yet, it is 
remarkable how well the scaling Ansatz eqs.(^,^3|) fits the numerical data. 



0.8 




2 4 6 8 10 12 



e(lo g (t)) 1/(1+a) 

FIG. 8. Distribution Go(e) as function of e-^/log(t) for case B (distribution with a maximum) and as function of elog(f) for 
case C (monotonically decreasing distribution). Times are t = 2 k with k = 6, 8, 10, 12, 14, 16. 
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0.5 



£ (lo 2 (0) 1/(1+a) 



FIG. 9. Distribution Gi(e) as function of ey^log(t) for case B (distribution with a maximum) and as function of elog(f) for 
case C (monotonically decreasing distribution). Times are t = 2 fe with k — 6,8, 10, 12, 14. 



VII. A METHOD TO DETERMINE THE THRESHOLD ENERGY SCALE e 



In this section we are interested in the following question. Is there a general method to determine the energy scale 
e* without having any precise information about the adiabatic modes present in the system ? In the previous sections 
we addressed this question by proposing an adiabatic scaling Ansatz to the dynamical equations. Here we propose a 
general method to determine the energy scale e* from first principles without the necessity of knowing the nature of 
the slow modes present in the system. Obviously for models such as the standard BG model this energy scale has no 
role since we know from the beginning that relaxation takes place on a single energy scale. 

Consider the following quantity P(AE) defined as the normalized probability density of having a first accepted 
energy change AE at time t. Let us consider the case of zero temperature where this probability density is defined only 
for AE < 0. If Q(AE) denotes the probability of proposing an energy change at time t (the move is not necessarily 
accepted), it is easy to show that P and Q are proportional to each other: 



where A — J_ Q(AE)dAE is the acceptance rate. The expression for Q(AE) (and therefore P(AE)) can be exactly 
computed. Note that computing Q(AE) yields all information about the statistics of energy changes, in particular the 
evolution equation for the energy Q On the contrary, given the time evolution for the energy this does not necessarily 
yield the distribution Q(AE). For the DBG this function can be exactly derived (its derivation is shown in the 
Appendix D). Here we quote the result, 



P(AE) 



Q{AE) 



6(-AE) 



(52) 



.4 



P(AE) 



J °° de 50 (e) 9l (e - AE) + (1 - P ) 9l (-AE) 



9{-AE) 



(53) 



.4 



with 




(54) 



Actually, in equilibrium at finite temperature Q(AE) satisfies detailed balance Q(AE) — Q(— AE) exp(— (3AE). 
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Using the scaling Ansatz eq.(|f)|) we obtain the simple scaling scaling relation, 

P(AE) = ^P(^) (55) 

A collapse of different P(AE) for different times can be used to determine the time evolution of e*. In figure ([h]) 
we show the scaling of P(AE) for the model B for N = 10 4 and different times t = 10 2 , 10 3 , 10 4 , 10 5 . Starting from a 
random initial configuration, statistics has been collected over approximately 30000 jumps for every time. In figure ( p"l| ) 
we check the scaling relation ( |55| ) plotting P{AE)e* as function of AE/e* where we have taken e* ~ l/-\/log(i). Note 
also that the range where P(AE) is finite corresponds to the region where e ~ e*. In figure ( [Tl| ) this corresponds to 
e* ~ 6/ ^/log(f) in agreement with what was observed in figures (||,^|) ■ 




1 r 




AE (log t) 1 



FIG. 11. Scaling plot for P(AE) / y/\og(t) versus AEy/\og(t) for different times t = 10 2 , 10 3 , 10 4 , 10 5 



The scaling works pretty well showing how this method could be used to guess the time evolution of the energy 
threshold e* in general glassy models in those cases where there exists an energy threshold. 



VIII. CONCLUSIONS 



In this paper we have considered a new solvable glass model, the disordered Backgammon model (DBG). The new 
ingredient of this model is that each box has associated a positive random energy which is obtained from a distribution 
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g(e). Again, similarly to its predecessor (the Backgammon model (BG)), the model displays slow relaxation due to 
the presence of entropic barriers. Actually, it turns out that the relaxation at T = of the number of empty boxes and 
all other occupation probabilities Pk is exactly the same as the original BG model, and, in particular, independent 
of the disorder distribution 17(e). In general, the relaxation of other quantities such as the energy and other disorder 
dependent observables, displays an asymptotic relaxation which depends on the statistical properties of 17(e) in the 
limit e — * 0. In the asymptotic long-time regime, relaxation takes place by diffusing particles among boxes with the 
smallest values of e. Therefore the asymptotic decay of the energy as well as that of other observables only depends 
on the exponent a which describes the limiting behavior g(e) — > e Q . The original BG model is recovered as a special 
case in the limit a — > 00. 

We have written integral equations for the densities <7fc(e). These equations form a hierarchy of dynamical equations 
which can be closed by introducing a suitable generating function. We focused on the solution of this hierarchy in 
the particular case of zero temperature. In this case the analytical solution of these equations proceeds in two steps. 
Firstly, the equations for the occupancies P& are exactly the same as in the original BG model and they can be solved 
by using known analytical methods. Secondly, this information is used to guess an adiabatic solution for the <7fc(e) in 
terms of a time-dependent energy threshold e* . 

All densities are out of equilibrium but admit a scaling solution of the type (^(e/e*), with the condition that, for 
e > e*, gk decays to zero fast enough to guarantee that dxcjk(x) = 1. That means that for energies above the 
threshold the modes are almost completely thermalized. 

This Ansatz solution yields two types of leading behaviors: on the one hand, they yield the asymptotic long-time 
behavior of e*,z* at zero temperature; on the other hand they produce the low-temperature behavior of the values 
of e*,z* in the stationary equilibrium limit t — > 00. The adiabatic approximation is nothing else but stating the 
validity of the complementary description of these two very different regimes. On the one hand, the equilibrium 
regime, where first the limit t — > 00 is taken and later T — > 0. On the other hand, the far from equilibrium regime, 
where the limit T — > is taken first and later t — > 00. The commutation of these two limits allows for the interchange 
between different variables such as energy, temperature and times when expressed in terms of their asymptotic leading 
behavior. Knowing the leading behavior of the quantities z* and e*, dimensional reasoning as presented in Appendix 
B yields the leading behavior of (3 e s which turns out to be proportional to z* /e*. The DBG model offers a scenario 
where there are two energy sectors separated by the energy scale e* which have very different physical properties. 
These two sectors manifest in the behavior of observables like the probability densities gk (e) where the time dependent 
threshold e* separates the equilibrated modes (e > e*) from the modes which stay off-equilibrium (e < e*). In the 
off-equilibrium regime entropic barriers are typically higher than the time dependent barrier at the threshold level e* . 
For e > e* barriers are lower and equilibrium is fastly achieved. 



IX. DISCUSSION 



What is the interest of the existence of this energy scale in the context of real glasses? One of the most interesting 
properties of the present model is that it introduces, in a very simple way, the concept of a threshold energy scale. 
To our knowledge such concept has never been discussed in any one among the plethora of mean-field glassy models 
studied during the last years. In those cases one studies relaxation of global quantities which get contribution of 
all possible energy scales involved into the problem. In principle, nothing is wrong with that since macroscopic 
observables are those quantities which are always measured in the laboratory. The problem arises when tackling 
issues related to the violation of the fluctuation-dissipation theorem and concepts such as the effective temperature 
and partial equilibration. Usually an effective temperature is defined in terms of the measured dissipation (response) 
and fluctuations (correlations) in the aging regime. This effective temperature is supposed to quantify the amount 
of energy tr an s fer .when the system is put in contact with a thermal bath and behaves in several aspects as a real 
temperatureL3r 6 H 34 l. There is a fundamental problem with this definition which is the following. Suppose one takes a 
piece of silica well below T g (for instance, at room temperature). The piece of silica is not in equilibrium (actually, 
it is always relaxing even if the relaxation rate can be extremely small and unobservable at room temperature) so 
one would be tempted to claim that its effective temperature (that describing equipartition among the set of non- 
thermalized degrees of freedom) is around T g well above the room temperature. Obviously, if we touch a piece of 
silica, then the hand plays the role of a thermal bath at the room temperature. Therefore, why we do not feel the 
effective temperature which may be hundreds of degrees above the room temperature? Note that the content of energy 
of the glass is still very large. Indeed, if the glass suddenly crystallized it would liberate all its latent heato. There 
are two possible explanations to explain this discrepancy. The first one was analyzed in the context of the oscillator 
model and assumes that the thermal conductivity is so small that heat transfer is negligible such a short time-scalea. 
The other explanation is that, when touching the glass, we are not touching the slow collective degrees of freedom 
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which still contain a lot of energy but the thermalized degrees of freedom. These two explanations are not totally 
exclusive. Assuming the existence of an energy threshold e* such that, above e* , all collective modes are thermalized 
at the bath temperature and below that threshold they are off-equilibrium, this offers an explanation about why when 
touching a piece of glass we feel it at the room temperature. Our hand only couples to the higher energy degrees 
of freedom (phonons) and not to the (much hotter) collective excitations. Yet, we cannot exclude that even if we 
coupled the bath to the hottest collective degrees of freedom then the conductivity would be extremely small and no 
heat transfer would be measured. From a different point of view, the two different explanations for the small amount 
of heat transfer established between a bath and the "hot" glassy system reduce into a single one: the existence of a 
threshold scale e* is consequence of the highly different orders of magnitude for the conductivities in the two energy 
sectors. Future studies of other glassy models will better clarify this issued. 

Finally, we have proposed a method to determine the threshold energy scale e* by computing the general probability 
distribution Q(AE). Preliminary investigations in other glassy models show that this distribution provides a general 
way to determine the threshold scale e*. Moreover it gives interesting information about fluctuations in the aging 
state although future work is still needed to understand better its full implications in our understanding of the aging 
regime. 

Further investigations in this model will address other issues such as the measurement of effective temperatures. 
For instance, it would be interesting to understand how the effective temperature, defined as the temperature of the 
thermal bath which does not produce a net thermal current when put in contact with the system, depends on the 
energy sector e probed by the bath. By coupling the bath with a selected set of modes of energy e we can understand 
whether there is a single effective temperature T e g for all modes or rather, there is an e dependent temperature. Note 
that local equilibrium in this model is only valid in the energy sector e > e* and it could well be that there does not 
exist a well defined effective temperature in the other sector e < e*. Nevertheless, the most natural possibility is that 
thermal fluctuations in the off-equilibrium sector e < e* are determined by the effective temperature (see eq.(|3^)) 
which determines the relaxation rate of the slow collective mode Pq. Future studies should enlighten this and other 
related questions. 
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Appendix A: Occupation probability density equations. 

In this appendix we derive the equations of motion for the occupation probability densities for box-energy between 
e and e + de. First we start from the densities of having zero particle in a box of energy e. 

In the next table we list the processes contributing to the evolution of the occupation probability density of boxes 
containing zero particle. In the left column we show the processes involved in terms of occupation numbers of the 
departure box and of the arrival box. In the right column we write the correspondent contribution of a given process 
to the variation of the occupation density, Ago(e). 



occupation 


contribution to Ago(e) 


nd — 1 n a = 


<W<Wo [S(e - e d ) - 6(e - e a )} [l + ( e -««.-«) - l) 9(e a - e d )] 


n a > 




n d > 1 n a = 


-(l-5 nd>1 )S nat0 e-^5(e-e a ) 



The particle for which a jump is proposed is chosen in box d with probability n d /N . The arrival box is chosen with 
uniform probability X/N. 

The total difference per particle in the probability density of empty boxes of energy e is then 
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N N 



Using eqs. 



p=0 a=0 

<5n,,,i(l - S n a ,o)S(e - e d ) - (1 - <5„ d ,i)6„ a ,oe -/3ea <5(e - e a )} 
and the following identities, 

N 



j7^2S nafi e(e a - e) 



-/9(e a -e) 



1 



a=0 
N 



d=0 
N 



de' g (e') 
'de' 9l (e') 



-0{e'-e) 



-^^2S riatQ 6{e- e a )e 



-0e a _ 



5o(e)e 



-0€ 



we get the equation of motion for go(e) (namely eq. (gj)): 
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9 5Q (e) = Hm A^) =5i(£) 



1/N 



Pi 1-e 



-/3e 



rf6'.g (e') (e-^ e '- £ ) - l) 
)+ r&' 3l ( £ ')(e^ (e - e ' ) -l) 



(56) 
(57) 

(58) 
(59) 
(60) 



We then consider the evolution of the probability density for boxes containing one particle. In the following table 
we list the processes contributing to the evolution of such occupation probability density. 



occupation 


contribution to Agi(e) 


n d = \ n a = {) 


o~n d ,iO~n a ,o [S(e - e d ) - S(e - c a )] 


1 + (e-W«.-«J _ i) - e,,) 


n a = l 


-S ni ,i8 na ,i [8{e - e d ) + «5(e - e a )J 


n a >l 


-<$n d ,i(l _ ^tio.i _ ^«o,o)'5(e - ed) 


rid = 2 n a = 


Sn d ,2S„ a ,o [S(e - e d ) + S(e - e„)] e~ pea 


n a = 1 


8 nd ^5 na ,i [S{e - e d ) - S(e - e a )J 


n a > 1 


S nd ,%{l - *n B ,l _ <W,o)<!>( e - £ rf) 


nd > 2 n a = 


(1 - <5n d ,2 - S ni ,i)8 nat o5(e - e a )e~P €a 


n a = 1 


_ (1 _ ^4,2 - <W,l) < Wl*( e - e a) 



Departure boxes are chosen with probability n d /N. Arrival box are chosen with uniform probability 1/N. 
Using again eqs. (|[^) and eqs. (|58]-|60|) we are able to derive the equation of motion for the probability density of 
boxes with one particles and energy equal to e: 



dgije) 
dt 



2.92 (e 
5o(e) 



1 + J™ <kgo(e)e-f> e - P j - gi{e) 2 + J™ de'g (e') (er^~^ - l) 

fc' 5l (e')(e- /3(e - e ' ) -1 



Pi 1-e 



For densities of boxes with k > 1 particle the scheme of the contributions is presented is the following table: 
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occupation 


contribution to Afffcfe) 




— h k 1 n - — - £"* I 

lb S /V f6(j rh J- 


^n d ,/i<5ri a ,fe-i<5(e - e ) 




n a = k 




rid 


— k n a = 






0<n a = /Kfe-l 






n Q = k — 1 








rc a = k 


-S nd: ko\ la M [b\e - e d ) + 5(e - 


e a )\ 




n a > k 


- S n d ,k (l - T,h=0 S n a ,h) <K e - 




rid 


= k+l n a = 


-8n d ,k+lO~na,a$( e ~ e d)e pea 




< n a < k — 1 


0~n a ,k+lO~n a ,h<>(e - Erf) 




n. Q = fe - 1 


5n d ,k+i()n a ,k-i [S{e - ed) + S(e - 


-<*)] 




n a = k 


S ni ,k+i&na,k — e<i)-8(e- 


OJ 




n a > k 


8n d ,k+i (l - ELo 5 n a ,ft) <5(e - 


-^) 


rid 


> k n a = k — 1 


(i - Eftti <Wft) S na< k-iS{e - 


-o 




n a = k 







Combining all the contributions we obtain for g k (e) equation ( |26| 



6>t 



(fc + l)p fe+ i(e) 1 + 



1 + fe + k 



deg (e)e- l3e - P 
deg (e)e-^ - P ( 



+ 5fc-i(e) 



Appendix B: Ansatz for the adiabatic approximation. 

In this appendix we show that the Ansatz solution ( f40| ) is asymptotically a solution of the equations (|33-35) at 
T = yielding the leading behavior e* (|4l|). We start by rewriting (Effl in the following way 



A 5fc (e) 



(61) 



where AP k = P k - 5 k ,o, Ag k (e) = g fc (e) - 8 k ,og(e), M^) = a; 5 fc (x) and J °° dx g k {x) = J Q °° dx h k (x)/x = 1. Here we 
will perform the analysis for the case A: = 0. The equations for k > can be done in a similar fashion. Substituting 
this expression into eq.(B3f) we get 



dg (e) dAP 1 



Ot 



dt e 



-ha - - 



AP 0ul .e.de*(t) APi 



(e*) 2 



e 

-APi 



/n(~) / rf £ ' 5 ( e ')-AP / de'-M^) 



(62) 



AP ° W e ^ 
e e* 



where /i W stands for the first derivative of ho (x) . Note that the scaling function ho does not depend on time, hence 
there is no term in that expression. Now introduce (|2^) in the first term of the l.h.s of ( |62| ) and multiply the 
whole equation by -^pr to obtain, 



d\oe:(f*) 

AP X ho(x) + x h' (x) g J = AP 



AP 



o Jo 



de'g(e>) 



dx'g {x) 



eg{e) 



APo 



+ h (x) 



AP 



dx'gi (x) 
(63) 



where g k {x) = (h k (x) / x) . From this equation we can guess the scaling behavior of all quantities in the asymptotic 
large-time limit e* — > 0. In the sector e < e* we use g(e) ~ e a obtaining J Q de'g(e') ~ e a+1 . Assuming all terms of the 
same order, we get for e ~ e* 



APo - ( e *) a+1 
d\og(e*) 



AP X ~- 



&t 



(64) 
(65) 
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Using the standard adiabatic results Pq = 1 — l/z*, Pi = l/(exp(z*)) we obtain the results (|4|^|). Note that 
the set of equations for ht are still impossible to solve. Only in certain regimes such as e C e* it may be possible to 
obtain results. There is a set of equations which couples the different hk- But this set of equations is time independent 
and should yield all the scaling functions g^ (x) once appropriate treatment is taking of the amplitude constant which 
fixes the leading behavior of e* . 

We also consider, as an example, the case in which the probability distribution of the quenched disorder becomes 
exponentially high at high values of e and zero for low values, namely we choose 

g(e)=exp(-^j . (66) 

For this choice J* de'g(e') ~ — eexp (— 4) — ^ T (0, -j), where the generalized Euler function T(0, x) goes to zero as 
x — » oo. In order to estimate e* from eq. ( |63| ) we notice now that for Pi eq. (|6^) is still valid, while for AP we obtain 

AP ~-e*expf-4) , (67) 



eventually yielding 



e*(t) ~ — . (68) 

log(log(i)) 



Appendix C: Approach to equilibrium of the occupation densities 5/t(e). 

We present the equations of motions for the occupation densities in the asymptotic regime. The values of the 



densities are expanded to first order around their equilibrium values: = g k + 5g 



dSg (e) 
dt 



n(e)il+ I de'gl\e') 
^o(e) 
9o^) 



-0(e'-e) _ 



e -0t + peq ( X _ e -f3^ + z 



de 



' s e V) ( 



e - Pt _ e -/3e') 



(1 



-f3e\ 



de Sgo(e') ( 

X) 

de'6 9l (e') 



de' S gi (e') 



dSgi(e) 2 

= -6g 2 (e) - 6 9l (e) 



,-0e _ 



2 + J de'g (e') 

e-P< + PP (1 - e-^) + z I de'g Q (e') (e 

Jo 

g c \e) |z 2 e-^ J™ de' 6g (e') (l - e~^') + z j~ de' Sg Q (e') (. 

poo pe 

(1 - e"*) / de' 6 9l (e') - / de' S 9l (e') e~^~^ - 1 



e ^ - e"^ 



A: = 1 



dSg k {e) 
dt 



k + 1 

^.9fe+i(e) <5fffc(e) I 1 + - ) + ^fc-i(e) 



CW 1 



yk+l 



-/3e 



1 - 



de' ^o(e') (1-e^) 



A; > 1 



(69) 



k = , 



(70) 



(71) 



In the above equations (3 is the inverse thermal bath temperature and z is the equilibrium fugacity at that temperature. 
As T goes to zero (j3 — > oo, z(/3) — » oo) the equations for the first order perturbation to equilibrium can be closed: 



dSg Q {e) 
dt 



S 9l (e)- / de' [S gi (e) g(e') + g(e) d gi (e')} 



k = , 



(72) 
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dS 9l (e) 

at 

dSg k (e) 

at 



/OO 
de' [S 9l (e) g(e') + g(e) S 9l (e')} 

= Sg k (e) + Sg k -x(e) 



k = 1 , 
k > 1 . 



(73) 
(74) 



In order to estimate the relaxation characteristic time to equilibrium at low temperature we can expand eqs. (|6£ 
First we introduce the asymptotic threshold energy e*(T) as the energy discriminating between thermalized and not 
thcrmalized collective modes at temperature T. If we define it through the relation e*(T) = Tz(T) and we use the 
relation (22) obtained by doing a low T expansion then we get, 



e*(T) = z Q T— 

where z is the coefficient of the leading term of z(T) at low T (see eq.(|22]): z(T) — zoT^fe. 
Then we expand eqs. (|69|), take e ~ e* and introduce the following adiabatic Ansatz, 

Sg k (e) = g k {e) - g c k q {e) = A ^ fc ^^ g fc (-^) . 



(75) 



(76) 



Note that this solution is equivalent to the Ansatz eq. ( |40| ) introduced for the asymptotic dynamics at zero 
temperature but with a static e*(T) now replacing the dynamical threshold. Now consider eq.(|69|) for 5go(e). Because 
SP k = J deSgk(e) it can be shown that the slowest mode corresponds to k = 0, i.e. Sgo(e) 3> 5g k {e) for k > 0. 
Therefore the second term in the right hand side of eq.(^) dominates the first and the second terms. Introducing 
( f76| ) into eq.([39|) we get that the relaxation time behaves like, 



T eq (e*) OC 



(77) 



For e >> e* the relaxation time is much smaller, since those are the modes with lower energy barriers. 
Appendix D: probability distribution of proposed energy updates. 



occupation 


contribution to E' — E 


probability 


rid = 1 n a = 


—£d + fa 


gi(ed) go{e a ) 


n a > 


-tel. 


siM [g{t a ) - go{£ a )\ 


n d > 1 n a = 




5o(e a )^ E P n p [g(e d ) - gi(e d )] 


n a > 





[g{£a) -9oUa)]jfJ2p n p[9^d)~9l(e d )} 



The probability distribution Q(AE) of proposed energy updates is the average of all possible changes, each computed 
with its probability: 



Q(AE) =5{E'-E- AE) 



(78) 



where AE is the proposed update, E is the energy of the system before the updating and E' the energy afterwards. 
This means 



OO C CO 



OO C CO 



Q(AE) = de de' 9l (e) g (e') S(A + e- e') + de de' 5l (e) [g(e') - g (e')} S(AE + e) 
Jo Jo Jo Jo 

r CO rOG -i 

+ j o de j o de' g (e')- ^ n p [g(e) - 9l (e)} 5(A e') 

+ j o de J a de' [ 9 {e') - go(e')]-Y, n P W) - ffi(e) }S(AE) 

pOO 

= / de 90 {e) 9l (e - AE) + (1 - P ) g\ (—AE) 9(-AE) - (1 - P x ) g {AE) 0(AE) 
Jae 

+(1-P )(1-Pi) S(AE) 



(79) 
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The term with S(AE) is the term responsible for diffusive motion of particles. Such a contribution does not actually 
give any contribution to the relaxation of the system and therefore we will not consider it from now on. 
The probability distribution of accepted energy changes is given by 



P(AE) = 

where W(f3AE) is the Metropolis function 

W(j3x) -- 

The normalization factor is 



Q{AE)W{f3AE) 
A 



e -0x 

1 



if x > 
if x < 



A= / de de 



- p{£ - £,) 9 o(e) 9i(e') + go(e') <?i(e)l + (1 - P )P 1 + (1 - Pi) / deg^e)e~^ 



As T — > the distribution P becomes 



P(AE) = 6{-AE) 



/ °° de 90 (e) 9l (e - AE) + (1 - P )gi(-AE) 



with A = / °° de' /~ de ff o(e') 9i(e) + (1 - P )Pv 

The normalization factor A is actually the acceptance rate of the Monte Carlo dynamics: 



A = 



dx W((3x) Q(x) 



as it was defined i nl. 

Using the same notation we can write the energy evolution as 



dE 
~dt 



dx x W{f5x) Q(x) 



+ de' de 9l (e') g (e) 



PiE- de 9l {e) e 
Jo 



(e - e 1 ) + (1 - I\) / de g (e)e-^e . 



(80) 



(81) 



(82) 



(83) 



(84) 



(85) 



The r.h.s. of this equation can be equivalcntly obtained following the procedure presented in appendix A. Indeed, by 
definition of energy density, is 



dE 
~dt 



de e 



%)(e) 
dt 



(86) 



Inserting eq. (p3) in eq. (pG) we get eq. (pi) back. 
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